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ABSTRACT 

Wc investigate the modal properties of inertial modes of rotating neutron stars with 
the core filled with neutron and proton superfluids, taking account of entrainment 
effects between the superfluids. In this paper, the entrainment effects are modeled by 
introducing a parameter rj so that no entrainment state is realized at 7y = 0. We find 
that inertial modes of rotating neutron stars with the superfiuid core are split into two 
families, which we call ordinary fluid inertial modes (i°-mode) and superfluid inertial 
modes (i^-mode). The two superfluids in the core counter-move for the z^'-modes. 
For the i°-modes, kq — liniii^Q uj/n is only weakly dependent on the entrainment 
parameter rj, where fl and to are the angular frequency of rotation and the oscillation 
frequency observed in the corotating frame of the star, respectively. For the i^'-modes, 
on the other hand, |ko| almost linearly increases as rj increases. Avoided crossings as 
functions of rj are therefore quite common between i"- and i^-modes. We find that 
some of the ?^-modes that are unstable against the gravitational radiation reaction 
at rj = become stable when 77 is larger than rjcrit, the value of which depends on 
the mode. Since the radiation driven instability associated with the current multipole 
radiation is quite weak for the inertial modes and the mutual friction damping in the 
superfluid core is strong, the instability caused by the inertial modes will be easily 
suppressed unless the entrainment parameter 77 is extremely small and the mutual 
friction damping is sufficiently weak. 
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1 INTRODUCTION 

Since Andersson (1998) and Friedman & Morsink (1998) found that the r-modes of a rotating neutron star are unstable 
against gravitational radiation reaction even for infinitesimally small rotation speeds, many papers on the r-mode instability 
have appeared. As a result of these investigations, the importance of the instability in astrophysics has been pointed out 
mainly in the following two ways, that is, the instability is strong enough to limits both the spin of newly born neutron stars 
and that of old accreting neutron stars, and neutron stars in which the r-mode oscillation is excited will be a potential source 
of gravitational waves that are detectable by the LIGO and other detectors (Lindblom, Owen, & Morsink 1998; Andersson, 
Kokkotas, & Schutz 1999; Owen et al. 1998; Bildsten 1998; Andersson, Kokkotas, & Stergioulas 1999). For a recent review, 
see, e.g., Andersson & Kokkotas (2001). 

It has long been suggested that neutrons in the inner crust and neutrons and protons in the core of neutron stars are 
in superfluid states when the interior temperatures cool down below T ~ lO^K (e.g., Shapiro & Teukolsky 1983). Since the 
interior temperature of neutron stars is believed to cool down quickly by emitting neutrinos (e.g., Baym & Pethick 1979), we 
may expect that many observable neutron stars, except newly born ones, are likely to have a core that contains superfluids. 
Although one of the important ingredients that affect the oscillations of neutron stars is superfluidity in the interior, it was 
rather recently that Epstein (1988) pointed out the importance of superfluidity for the oscillations, suggesting the existence 
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of a new family of acoustic modes associated with superfluidity by means of the local dispersion relations derived for the 
case of non-rotating neutron stars. Later on, Lindblom & Mendell (1994) studied the /-modes of rotating neutron stars that 
have a core filled with neutron and proton suporfluids, taking account of cntrainmcnt effects between the supcrfluids, and 
showed that the ordinary fluid-like /- modes are not strongly affected by the superfluidity in the core. In the same paper, 
Lindblom & Mendell (1994) suggested the existence of superfluid-like acoustic modes, which they called so-modes, but they 
could not find the so-modes numerically. First numerical examples of the superfiuid-liko acoustic modes of neutron stars with 
the superfluid core were obtained by Lee (1995), who also showed that g-modes are not propagating in the core filled with 
superfluids. Andersson & Comer (2001) discussed the dynamics of superfluid neutron star cores examining local dispersion 
relations in detail, and they confirmed that g-modes are excluded from the superfluid core (Lee 1995). Applying their argument 
to the r-modes, they suggested that the r-modes of a superfluid neutron star core are split into two distinct families, that 
is, ordinary fluid-like r-modes and superfluid-like r-modes. Most recently, Prix & Rieutord (2002) have computed in great 
detail acoustic modes of non-rotating neutron stars composed of neutron and proton superfluids for a wide range of the 
entrainment parameter. Although their results are basically consistent with the results obtained by Lee (1995) and Andersson 
& Comer (2001), they have also shown that the assumption that \dvp — dvn\ ~ for the ordinary fluid-like acoustic modes and 
\pp5vp + pn5v„\ ~ for the superfluid-like acoustic modes is not necessarily helpful to distinguish the two families, where 5vp 
and 5vn are the Eulerian velocity perturbations of the proton and neutron superfluids, the mass densities of which are denoted 
as pp and p„. We also note that Sedrakian & Wasserman (2000) examined rotational effects on the fundamental modes in 
uniform density superfluid stars, applying the so-called tensor virial technique to normal modes. As for rclativistic non-rotating 
superfluid stars. Comer, Langlois, & Lin (1999) calculated superfluid-like acoustic modes and showed that the basic features of 
the superfluid modes axe quiet similar to those of the corresponding Newtonian superfluid modes. More recently, Andersson, 
Comer & Langlois (2002) have found avoided crossings between the superfluid- and the ordinary fluid-like acoustic modes in 
a relativistic superfluid neutron star as the entrainment is varied. 

To give an answer to the question whether the r-mode instability will survive the damping due to mutual friction in the 
superfluid core of cold neutron stars, Lindblom & Mendell (2000) computed ordinary fluid-like r-modes, taking account of 
the entrainment effects in the core. They showed that the mutual friction damping in the core is not effective to suppress 
the instability except for a few domains of the entrainment parameter rj. Recently, Lee & Yoshida (2002) have studied the 
r-modes in neutron stars with the stipcrfluid core with a different method, confirming most of the results by Lindblom & 
Mendell (2000). Lee & Yoshida (2002) have also found that, as suggested by Andersson & Comer (2001), the r-modes are split 
into ordinary fluid-like r-modes and superfluid-like r-modes, and that the instability caused by the superfluid-like r-modes is 
extremely weak and easily damped by dissipation processes in the interior. 

The purpose of this paper is to study modal properties of inertial modes in rotating neutron stars containing superfluids 
in the core. As shown by Lockitch & Friedman (1999), the existence of inertial modes in the lowest order of angular frequency 
Q of rotation depends on the structure of the so-called zero-frequency subspace of eigensolutions to linearized hydrodynamic 
equations. Naively speaking, we may say that, for inertial modes to exist in the limit of O — > 0, g-modes of a non-rotating 
star with Q = need to be degenerate to the zero-frequency subspace. Investigating the zero-frequency subspace of pulsation 
equations for the superfluid core, Andersson & Comer (2001) showed the existence of the zero-frequency g-modes, and hence 
the existence of inertial modes (see, also. Comer 2002). In this paper, we calculate global inertial mode oscillations in neutron 
stars with a core filled with neutron and proton superfluids, taking account of the entrainment effects between the two 
supcrfluids. For the dynamics of supcrfluids in the core, we employ a two-fluid formalism similar to that given in Lindblom & 
Mendel (1994). To obtain the inertial modes in the limit of ^ 0, we generalize the method devised by Lockitch & Friedman 
(1999) to stars composed of a superfluid core and a normal fiuid envelope. In §2 we present the basic equations employed in 
this paper for the dynamics of superfluids in the core, and in §3 dissipation processes considered in this paper are described. 
§4 gives numerical results, and §5 and §6 are for discussions and conclusions. 



2 OSCILLATION EQUATIONS 

We derive basic equations governing superfluid motions in the neutron star core in the Newtonian dynamics, assuming tmiform 
rotation of the star. The core is assumed to be filled with neutron and proton superfiuids and a normal fluid of electron. We 
also assume perfect charge neutrality between the protons and electrons because the plasma frequency is much higher than 
the oscillation frequencies considered in this paper (see, e.g., Mendell 1991a). Since the transition temperatures Tc ~ 10^ K to 
neutron and proton superfiuids are much higher than the interior temperatures of old neutron stars (see, e.g., Epstein 1988), 
we assume that all the neutrons and protons in the core are in superfluid states and the excited normal fluid components of 
the fluids can be ignored. 

In this study, we assume that the neutron and proton superfluids and the electron normal fluid in the equilibrium state 
are in the same rotational motion with the angular velocity f2 around the axis of rotation. The velocity of the fluid in the 
equilibrium state is therefore given by = iiip'^, where i^" is the rotational Killing vector. In a perturbed state, however, the 
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neutron and the proton superfluids can move differently from each other in the core, obeying their own governing equations. 
In order to describe the perturbations in the superfluid core, thus, we can choose the perturbed neutron density Spn, the 
perturbed proton density Spp, the perturbed velocity of the neutron Sv"^, and the perturbed velocity of the proton Svp as 
dynamical variables. Since the equilibrium star is axisymmetric about the rotation axis, we assume that the time and azimuthal 
dependence of the perturbations is given by exp{iat + imtp) with a being the oscillation frequency observed in an inertial 
frame, and m is an integer representing the azimuthal wave-number. The basic perturbation equations employed in this paper 
for the neutron and proton superfluids in rotating neutron stars are essentially the same as those given in Mendell (1991a) 
and Lindblom & Mendell (1994) and are given as follows: 

iuSpn + Va (pn(5lC) = 0, (l) 
lLu5pp + Va {Pv5Vp) = 0, (2) 

iiuSvl + 20. ^w^ViVs" + ^"{Sp^ + 5$) -20.^ {5vl - 5v'^,) VVt = , (3) 

Pn 

iLu ( Sv'i, H Sv'^ ) + 2Q I Svt H Sv'^ ) Vtv?" 

-ip J \ J 



+\/''iSpp + ^Sp,+ [l + ^]5^]-2n^ {5v'„ - 5v;) = , (4) 

\ TTlp \ ITT-P J J pn 

where lu denotes the corotating frequency, defined by tj = a + mQ, Va means the covariant derivative in fiat three dimensional 
space, and SQ stands for the Eulerian change in the physical quantity Q. Ifere, mp and me are the proton and the electron 
masses. In the mass conservation equations Q and l(2j, pnSV^ and ppSVp mean the mass current vector of the neutron and 
of the proton, respectively, and can be expanded in terms of the perturbations of the superfiuid velocities as follows: 

PnSV^ — pnn Sv^ + pnp 5Vp , (5) 

PpSVp = ppn Sv"^, + ppp Svl , (6) 

where the relationships p„ = pnn + pnp, pp ~ ppn + ppp, and pnp = Ppn must be satisfied due to the Galilean invariance of 
equations (0 and The perturbed superfiuid velocities are then given in terms of SV,f and SVp as 

<s< = ^ sv: + ^ 5v: , (7) 

Pn Pn 



5y;=P2L5V: + ^SVp\ (8) 
Pp pp 

where 



2 2 
PppPn Pnnpp pnpPnPp 
Pll = -2 . P22 = r^— , Pl2 = P21 = , (9) 

and 

P — pnn ppp Pnp ppn- (10) 

Note that pn + pi2 = pn and p22 + P21 = pp- In Euler equations JHJi and Q, Sv"^ and (S"l> are the perturbed velocity of 
the electron and the perturbed gravitational potential, respectively, and 5p,n, Spp, and Spe stand for the perturbed chemical 
potentials per unit mass for the neutron, proton and electron, respectively. Because of the assumption of the perfect charge 
neutrality in the proton and electron plasma, we may have 

SVp — Sv^, and Spp/rUp = Spe/rrie, (11) 

where pe is the electron density. The perturbed gravitational potential is determined by the Poisson equation 

V^(5<E> = iivGSp, (12) 

where p — pn + Pp + Pe and G is the gravitational constant. 

Using a variant of Gibbs-Duhem relation, the pressure perturbation is given by 

5p = p„5pn + PpSp , (13) 

where Sp = 5pp + ^ 5pe- Here, we have ignored the entropy carried by the electron normal fluid. This treatment is justified 
for a star whose temperature is low enough. Note also that the superfluids carry no entropy. 

To obtain a relation between the densities p„ and pp and the chemical potentials p„ and p, we begin with writing the 
energy density e as 
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e = e{pn,Pp) , 

with which the chemical potentials are defined as 

(Pn,Pp) = {de/dpn)p^ , Pp ipn,Pp) = {de/dpp)^^ 
If the chemical potential of the electron is given by 



we have, due to charge neutrality condition pe 
5/in \ _ "Pu 'Prz \ ( Spn 



ppirie/mp, 



\ P21 P22 J \ Spp 



where 

Pu = 



dpn 



P12 = 



Pp 



9p, 



9/i 
dpp 



We write the inverse of equation l|17|l as 



Spn 
Spp 

where 



Qii 
Q21 



Qii 
Q21 



Q12 
Q22 



Q12 
Q22 



Pii Pi 
P21 V2 



The curl of the Euler equations, which is the vorticity equation, is useful for perturbation analysis of rotating 
is given in terms of SV^ and SVp by 



abcT-7 



lUJ — Qcd + 2^2 Vdiy3c I + io;— Qcd&Vp 

Pn I pn 



= 0, 



(14) 
(15) 
(16) 

(17) 
(18) 

(19) 

(20) 
stars, and 

(21) 



abcr—! 



P22 ^ rUc 
Pp rrip 



20, [ I 



Vd'i^c ] SVp + iuj— QcdSV^ 

' pp 



0, 



(22) 



where gat and mean the metric tensor and the Levi-Civita tensor in the flat space. 

Although we have made no assumptions for rotation velocity up to now, let us restrict our consideration to slowly rotating 
stars in order to obtain inertial modes to the lowest order in the stellar rotation frequency fi. In the limit of slow rotation, 
then, we assume perturbations to obey the following ordering in Q: 



<S1C = 0(1), SVp = 0(1), Lj/n = 0(1), Sp„ = 0(Q.), Spp = 0{n), and S$ = 0{Q.) . 



(23) 



Examining the zero-frequency subspace of the perturbation in a non-rotating superfluid core, Andersson & Comer (2001) 
showed that solutions obeying the ordering 1231 may be allowed in the equations describing hydrodynamics in the superfluid 
core. The existence of solutions obeying 12311 is similar to that in a barotropic normal fluid star (Lockitch & Friedman 1999). 
Those solutions can be interpreted as infinitely degenerate g-modes of zero- frequency if we consider no rotational effect. It is 
convenient to introduce dimensionless frequency kq in the limit of S7 — > 0, which is defined by 



Ko ~ lim Lu 10, . 

Then, basic equations describing the lowest order inertial modes are given by 
v. (p„<50 = 0, V. (pp<5V7) = 0, 



(24) 



(25) 



ahcr—! 



Pn n- -lv7 \ jTT/d 1 Pl2 s:trd 

Qcd - 2i«:o Vd((9c oV„ H QcdoVp 

pn I Pn 



0, 



(26) 



a abc-r-7 



P22 _|_ ma 



- 2tK^' ( 1 + ^ Vd^c SVp^ + ^ g^dSV:^ 
m„ J J pp 



= 0. 



(27) 



Pp iit^p I \ '"'p 

Note that q'^ (qp) has only two independent components because VaQn ~ (Vagp=0) is automatically satisfied. 

/From similar consideration to that for the superfluid core, basic equations for pulsations in the normal fluid envelope 
are given by 
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Va ipSv") = , 



(28) 



= e'''"'\/t [(ged - 2iKo'\/a<pc) Sv"] = , (29) 

where Sv"' denotes the perturbation of the fluid velocity. 

The perturbations in a uniformly rotating star are expanded in terms of spherical harmonic functions with different Vs 
for a given m (e.g., Lee & Saio 1986). The vector 5Vx is therefore given by 



i>|m| 



E 



^.^w|vr(.,^)-.r,f(.)-i^^y-(.,^) 



(30) 
(31) 



!>jm 



(r) (e, V') + (r) sin {9, ^) 



(32) 



where I — \m\ + 2{k — 1) and I' = I + 1 for even modes, and I = \m\ + 2A: — 1 and I' — I — 1 for odd modes, where 
fc = 1, 2, 3, Here, we have employed the constituent index X which can be either n, p or o, where o means the 

velocity field in the normal fluid envelope, and the normal spherical polar coordinate (r, 9, ip) is used. In this paper, the even 
and odd modes are respectively characterized by symmetry and antisymmetry of their eigenfunctions with respect to the 
equatorial plan. Substituting these expansions into equations 1)25^ through l|27|l . we obtain oscillation equations given as a set 
of simultaneous linear ordinary differential equations of the expansion coefficients (see Appendix), which is to be integrated 
in the superfluid core. The oscillation equations solved in the normal fluid envelope are similar to those given in Lockitch & 
Friedman (1999). 

To obtain a unique solution of an oscillation mode, solutions in the superfluid core and in the normal fluid envelope are 
matched at the interface between the two domains by imposing jump conditions given by 



5v: = 5v; = sv: 



(33) 



The boundary conditions at the stellar center is the regularity condition of the perturbations 5V^ and 5Vp. The boundary 
condition at the stellar surface is Ap — 0, where Ap is the Lagrangian change in the pressure. 

For numerical computation, oscillation equations of a finite dimension are obtained by disregarding the terms with I 
larger than Imax ~ \m\ + 2kmax — 1 in the expansions of perturbations such as given by (1301 to 1321 1. where Imax is determined 
so that the eigenfrequency and the eigenfunctions are well converged as Imax increases. We solve the oscillation equations of 
a finite dimension as an eigenvalue problem with the scaled oscillation frequency kq using a Henyey type relaxation method 
(see, e.g., Unno et al. 1989). 



3 DISSIPATIONS 

In the previous section, we consider no dissipation effects. Thus, the conserved energy E of oscillations exits and is given by 
(MendeU 1991b) 

7^2 



+ I f (pSvlSv'' + ( ^ + 5$ ) 5p' ] (34) 



where g is the determinant of the metric, and asterisk (*) denotes the complex conjugate of the quantity. Here, dPx and 
d^x means the integration over the superffuid core and the normal fluid envelope, respectively. In this equation, we have 
employed new variables, deflned by 

Sv"- = (pn(5< + pj,5Vp)/p= {pnSVn + ppSVp) / p , 

Sw-^ = Svl - 5< = ^ {SV: - 5V^) , (35) 

Spi — Spn , Sp2 = 5pp . 

Note that pSv'^ is interpreted as the total momentum current vector of perturbation. In reality, however, there are a lot 
of dissipation mechanisms in real neutron stars. Thus, the oscillation energy E is not conserved in real neutron stars. In 
other words, the amplitude of an oscillation mode should damp or grow depending on considered dissipation mechanisms. 



6 S. Yoshida and U. Lee 



The stability of an oscillation mode of a star is then determined by summing up all contributions from various dissipation 
mechanisms. In this study, we consider the contributions from gravitational radiation reaction, viscous processes, and mutual 
friction in the superfluid core. Then, the averaged energy loss (or gain) rate dE/dt of a normal mode in a rotating neutron 
star may be given by 

^= -cruj ^ Nia^' \SJi,,j,\'^ -2 / ■qsScj''^* 5cjab^/g d^x 

l'>\m\ JVs + Vn 

— 2Vl / p„B„ i — — I (gab — ZaZt)5w''*Sw^y^ d^x 

Jvs \PnPpJ 

- i:(f)„,„-(f),-(f)„' 

V>\m\ 

where z°' is the unit vector that parallels the rotation axis, defined by Za = eai,cV'i/j'^/2. 

The terms {dE/dt)cj,i on the right hand side of equation 13611 denote the energy loss (or gain) rates due to gravitational 
radiation reaction associated with the mass current multipole moment 5 Jim, where 

&Jim = , / r^^ pSvae'''"'\7 bYr'*\/ cry^ d^ X (37) 



Vs+v„ 



and 

_ 47rG (/ + !)(/ + 2) 

' c2'+iZ(/-l)[(2Z + l)!!]2' ^ ' 

and c is the velocity of light (Thorne 1980, Lindblom et al. 1998). 

The term {dE/dt)s is the energy dissipation rate due to shear viscosity, and rjs is the shear viscosity coefficient, and Saab 
is the traceless rate of strain tensor for the perturbed velocity field, which is defined by 

SUab = ^ ^Va^tb + VbUa — ^Qab^ cU"^ , (39) 

for the velocity field ii". In this paper, we ignore the contribution from the bulk viscosity, which is important only for newly 
born hot neutron stars without superfluids in the core. The shear viscosity coefficient we use in the superfluid core is 

Vs = 6.0 X 10- (j-^) ' ('4^) ' g/cm s (40) 



^10i5g/cm3 J \ T ^ 

(Cutler & Lindblom 1987, Flowers & Itoh 1979), and that in the normal fluid envelope is given by 



- loW ) ["^ ) ^/^"^ ^ ^''^ 

(Cutler & Lindblom 1987, Flowers & Itoh 1979). The rate of strain Saab is evaluated by using Sv^ in the superfluid core, and 
in the normal fluid envelope (Lindblom & Mendell 2000). 

The term [dE/dt)MF is the energy loss rate due to the mutual friction in the superfluid core, and the dimensionless 
coefficient _B„ is given by (Mendell 1991b) 

= 0.011 ^P^(eZL\'' (PlA' ( ) . (42) 



9/4 



pn\Pv J \Ppv J \10^*g/cm 

The mutual friction is a dissipation mechanism inherent to a rotating system of superfluids, and it is caused by scattering of 
normal fluid particles off the vortices in the superfluids. Since we have assumed perfect charge neutrality between the electrons 
and protons, we consider scattering between the normal electrons and vortices of the neutron superfluid (Mendell 1991b, see 
also Alpar, Langer, & Sauls, 1984). 

As is indicated by the first term on the right hand side of equation 13611 . if the normal mode has an oscillation frequency 
that satisfies —au > 0, the oscillation energy E in the corotating frame, in the absence of other damping mechanisms, increases 
as a result of gravitational wave radiation, indicating instability of the mode (Friedman & Schutz 1978). 

The damping (or growth) time-scale r of a normal mode may be given by 

1 1 (dE\ 1 11 

!'>|mj 

where = —2E/ (dE/dt)i. For the inertial modes, it is convenient to derive an extrapolation formula of the time-scale r given 
as a function of SI and the interior temperature T (e.g., Lindblom et al. 1998, Lindblom & Mendell 2000): 
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i'flll ^G,/..' U^p j +r«l^ T ) +r^^^[nGp) ' ^ > 

where p = M/(47rilV3)- Here, M and are the mass and radius of the star. The quantities Tqj i,, t^, and assumed 
to be independent on Q, and T. 



4 NUMERICAL RESULTS 

Following Lindblom & Mendell (2000), we employ a polytropic model of index N = 1 with the mass M = IAMq and the 

radius R — 12.57km as a background model for modal analysis. The model is divided into a superfluid core and a normal fluid 
envelope, the interface of which is set at p = ps = 2.8 x 10^'*g/cm^. In the normal fluid envelope, we use the polytropic equation 
of state given by p = Kp^ both for the equilibrium structure and for the oscillation equations. The core is assumed to be filled 
with neutron and proton supcrfluids and a normal fluid of electron, for which an equation of state, labeled A18+5i;+UlX 
(Akmal, Pandharipande, & Ravenhall 1998), is used to give the thermodynamical quantities used for the oscillation equations. 
For the mass density coefficients p„„, Ppp, and p„p in the core, we employ an empirical relation given by (see Lindblom & 
Mendell 2000) 

Pp/p w 0.031 + 8.8 X 10""p, (45) 
and a formula given by 

Pnp = -VPn, (46) 

where r; is a parameter of order of ~ 0.04 (Borumand, ,]oynt, & Kluzniak 1996). In this paper we call rj the ontrainment 
parameter. In this study, only the modes having m = 2 are examined because they are the most significant for the gravitational 
radiation driven instability. For the range of the entrainment parameter rj, we choose —0.01 < rj < 0.04. Although negative 
values of rj might not be suitable for physical models of neutron stars, we include in our calculations negative r] to demonstrate 
clearly the rj dependence of the inertial modes around 77 ~ 0. 

We find that inertial modes in the superfluid core are split into two families, which we call ordinary fluid inertial modes 
(i°-modc), and superfluid inertial modes (i^-modc). This mode splitting is common in other oscillation modes in the superfluid 
core, because the dynamical degree of freedom in the superfluid core simply doubles (see, e.g., Epstein 1988; Lindblom & 
Mendell 1994; Lee 1995; Andersson & Comer 2001; Prix & Rieutord 2002). As a subclass of inertial modes, the r-modes 
of neutron stars with a superfluid core are also split into ordinary fluid-like r-modos and superfluid-like r-modes, which we 
respectively call r°-modes and r^-modes as in Lee & Yoshida (2002). Since basic properties of the eigenfunctions (5v" (or Sw'^) 
of an inertial mode are quite similar to the eigenfunctions (5v" of the corresponding inertial mode in the normal fluid star 
without core superfluidity, we employ the same classification scheme as that employed by Lockitch & Friedman (1999) and 
Yoshida & Lee (2000a) for inertial modes in the normal fluid star. In this scheme, inertial modes are classified in terms of the 
two quantum numbers m and lo — \m\, and the scaled frequency /to = limn->ow/f2- For given m and lo — \m\, the number 
of inertial modes of one kind is Iq — \m\. Even parity inertial modes have oven Iq — |m|, and odd parity inertial modes have 
odd lo — \m\. Note that the r-modes are classified as inertial modes with lo — \m\ = 1 (e.g., Yoshida & Lee 2000a). It is 
interesting to note that, although there exist three quantum numbers to classify the normal modes of non-rotating stars, a 
third quantum number may be missing for inertial modes. Usually, the ratio Es/E of the oscillation energy Es in the core to 
the total oscillation energy i5 is a good indicator for a distinction between i"- and i°-modes, since we find Es/E <0.5 for the 
i°-modes and Es/E >0.9 for the i°-modes for the models composed of a superfluid core and a normal fluid envelope. 

Lot us first consider the case of 77 = 0, that is, the case of no entrainment. In Table 1, kq's for i°- and i^-modes of 
low /o — are tabulated together with no's for the corresponding inertial modes in the iV = 1 polytrope model without 
superfluidity in the core, which we simply call i-mode (see, Lockitch & Friedman 1999; Yoshida & Lee 2000a). Our values of 
Ko for the i°- and i^-modes are in good agreement with those computed by Lee & Yoshida (2002) with a different method. 
Note that we have shown five digits for numerical values of ko to indicate that our numerical results include at most a few 
percent numerical error. We can see that ko's of paired i°- and i*-modes have quite similar values to that of the corresponding 
inertial mode in the N — 1 polytropic normal fluid model. In order to demonstrate model properties of the inertial modes, 
we exhibit the expansion coefficients Sf , T^i , SJ" , and T^f of the eigenfunctions Sv"^ and Jui", defined by 

^{p„Sr + PpSf)/p, Tr,={p^T{} + ppTf,)/p, (47) 

ST = P^Pp{ST^S^)/p\ Tlf ^p^pp{TP-Tf,)/f. (48) 

Note that Sf and T/f vanish in the normal fluid envelope because neutrons and protons co-move. In Figures 1 through 4, 
the first three expansion coefficients Sf, TjT, 5f, and are shown for the i°- and i*-modes of lo — \m\ = 2 which have 
Ko = —0.5685 and ko = —0.5115, respectively. Here, the amplitude normalization given by T^+i = 1 has been used. It is 
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noted that Sv"^ of the i°-modes and Sw" of the i^-modes resemble Sv"^ of the corresponding inertial modes of the normal fluid 
star model (compare the figures with Figure 1 of Yoshida & Lee 2000a). It is also found that although the amplitudes of 

and 5w'^ arc similar to each other for the i°-modcs, the amplitude of Sw'^ dominates that of Sv" for the i^-modes, which 
indicates that the two superfluids in the core basically counter-move for the i^-modes. 

In Figures 5 to 8, ko's of inertial modes of low lo — \m\ are plotted as functions of the entrainment parameter rj, where 
Figure 5 is for even parity modes with positive ko, Figure 6 for even parity modes with negative kq, Figure 7 for odd parity 
modes with positive kq, and Figure 8 for odd parity modes with negative kq. In Figure 7, we have also given kq of the r^-mode, 
the locus of which as a function of rj is in good agreement with the locus calculated by Lee & Yoshida (2002) with a different 
method. Note that we cannot obtain r°-modes by using our numerical method, the reason for which is not clear (sec the next 
section) ^. From Figures 5 through 8, we find that, although kq for the i°-modes is only weakly dependent on the parameter 
rj, \ko \ for the i^-modes increases almost linearly as rj increases. Because of these different dependences of kq on the parameter 
rj, frequent avoided crossings are inevitable between i°-modes and i*-modes. As shown by Figure 7, for example, the i^-mode 
of lo — \m\ = 3 which has ko = 1.3865 at 77 = suffers three avoided crossings with different i°-modes as rj is increased from 
—0.01 to 0.04. Since kq's of paired i°- and i^-modes at 77 = are close to each other, they always make an avoided crossing 
around rj = 0, which is also shown by Figures 5 through 8. 

In Figures 9 through 12, we display the first three coefficients Si , T^,, SI" , and for the i°- and i^-modes of lo — \m\ = 2, 
which respectively have ko = —0.5403 and ko = —0.6103 at 77 = 0.01. Here, the amplitude normalization given by T^+i = 1 
has been used. Note that these two inertial modes have kq = —0.5685 and ko — —0.5115 at 77 = 0. Comparing these figures 
with Figures 1 to 4, we observe that the basic properties of the eigenfunctions of the inertial modes do not depend on the 
entrainment parameter rj, except at avoided crossings. 

In Tabic 2, wo tabulate ko and the scaled damping (or growth) timescales rf for the i°- and i^-modcs with positive kq for 
77 = 0, 0.02, and 0.04. We find that the growth timescales Tqji, for all the i°- and i°-modes in the table are longer than that 
of the r-mode of the normal fluid neutron star, which means that the instability of the i°- and z*-modes caused by the current 
multipole radiation is not strong. This result is the same as that for the inertial modes in a normal fluid star. It is interesting 
to note that the two i^-modes in the table having ko = 2.1661 and ko = 2.4083 at 77 = 0.04 are stable against gravitational 
radiation reaction since they have positive Tqj ;/ . The reason for this can be understood as follows. As mentioned in the last 
section, the gravitational radiation driven instability sets in when the frequency of the mode satisfies a(a — 7nf2) < 0, which 
leads to the instability criterion given by < /to < 77X for positive m. For inertial modes of normal fluid stars, it is known that 
all inertial modes having positive ko for 771 > satisfy the instability criterion (Lindblom & Ipser 1999; Lockitch & Friedman 
1999; Yoshida & Lee 2000a). For inertial modes in the superfluid core, on the other hand, since the frequency |Koj of the 
i^-modes monotonically increases as the value of 77 is increased, the ko of some positive-frequency i^-modes exceeds the limit 
m when 77 gets larger than a critical value rjcrit, the value of which depends on the mode (see, Figures 5 and 7). Thus, some 
of the positive-frequency i''-modes that are unstable satisfying the instability criterion at r; = become stable invalidating 
the criterion for 77 larger than rjcru- As for the effects of mutual friction, we observe that the damping timescales Tmf for the 
inertial modes are much shorter than those due to other dissipation processes unless mutual friction between the superfluids 
is extremely weak, that is, the parameter 77 is extremely small. Note that the amount of dissipation caused by the mutual 
friction is roughly proportional to Bn\Sw'^\^ where -B„ oc 77^. Accordingly, if we assume 77 ~ 0.04 as a typical value for the 
parameter, the dissipation due to mutual friction in the superfluid core is very effective to damp out the i°- and i^-modes. 



5 DISCUSSION 

As mentioned in the previous section, with our method of calculation we fail in obtaining ordinary fluid r-modes, which are 
called ^"-modes in Lee & Yoshida (2002). However, this docs not necessarily mean that our basic equations (A2) (A7) do not 
contain such r-mode solutions. To prove this, let us consider purely axial parity solutions, whose angular behavior is given by 
a purely axial harmonic function having I' = m. In this section, we also assume for simplicity that the densities p„, pp, and 
pi2 are constant and do not depend on r. In virtue of the assumption of purely axial mode, we can reduce our basic equations 
to the following equations (see. Appendix): 
For the superfluid core, 

^ Recently, wc have generalized our formalism to the case of general relativistic superfluid neutron stars. In that study, our numerical 
code for Newtonian dynamics have been straightforwardly extended into the one for general relativistic dynamics. With this relativistic 
version of the code, both r°- and r'-modes can be obtained even in almost the Newtonian regime with small values of the relativistic 
parameter M/R. Fundamental differences in the two formalisms reside only in the master equations. As shown in equations (49)-(52), 
the basic equations for neutron and proton fluids in the Newtonian case get to degenerate in lowest order of i1 when r°-modcs are 
considered, while the basic equations in the relativistic case don't even if we are concerned with r°-modes (see, Yoshida & Lee 2003). 
We therefore guess that the reason why r°-modes cannot be obtained in the Newtonian case is this possible degeneracy of the two sets 
of superfluid equations in Newtonian dynamics. 
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um{m + 2) J"+i {'■^^™ 
um{m + 2) J™+i 



1 - 



rn{m + 1) 



Pl2 
/0l2 

m(m + 1) y pp 
For the normal fluid envelop 
d 



um{m + 2)J™+i |r 



m{m + 1) 



r° = 



m-l)T^}=0, 
m- 1)T^| =0, 

-T") = 0, 



T'' 1 = 



=0, 



(49) 
(50) 
(51) 

(52) 

(53) 
(54) 



where i/ — 2Q/a and J™ = [{l'^ — m^)/(4P — l)]^'''^. Here, we have ignored the terms due to m^/mp for simplicity. The 
condition for the continuity of the total momentum flux at the interface between the superfluid core and the normal fluid 
envelope is given by 



= {pnT^ + PpT^)/p a.t r = rs 



(55) 



where Vs is the radius of the superfluid core. Because the factor vm{m + 2) J^+i in equations H49^ . H5()|l and l|5H|l does not 
vanish for v and m > 2, we can obtain solutions to differential equations 1491 . (I50I I. and 15311 as 



T° = B{r/rs 



(56) 



where An, Ap, and B are constants, which are to be determined by boundary condition (55) and constraint equations (51), 
(52), (54). Substituting the solutions (56) into algebraic equations (51), (52), (54) and (55), we obtain 



f-pl2/pn pl2/pn 

pw/pp f-pu/pp 



An 

A.r> 



= W 



= 0, 



fB = 0, 
and 

B = {pn{rs)An + Pp{rs)Ap)/p{rs 



(57) 
(58) 

(59) 



where f — 1 — u/{m +1). For non-trivial solutions to exist in the core, we must have detM^ — 0, which leads to / = or 
/ = pi2{pn + Pp)/{pnpp)- For the case of / = (i.e. kq = 2/(m + 1)), solutions may be given by 



An — Ap 



B. 



(60) 



This solution obviously corresponds to the r-°-mode since the oscillation frequency ko is exactly the same as that of the 
r°-mode, and the two superfluids flow together in the core. For the case of / = pi2{pn + pp)/{pnpp) (i.e., kq / 2/(m + 1)), on 
the other hand, solutions are given by 



(p„(rs)A„ + pp{rs)Ap)/ pivs 



B = 0. 



(61) 



This solution corresponds to the r^-modes found by Lee & Yoshida (2002), for which the two superfluids counter- move in the 
core (see, also, Andersson & Comer 2001). It is to be noted that the solution (60) for arbitrary rj and the solution (61) for 
77 = can be derived even without the assumption that the densities pn, Pp, and pi2 do not depend on the radial distance r. 



6 CONCLUSION 

We have discussed the modal properties of inertial modes of rotating neutron stars with the core filled with neutron and 
proton superfiuids, taking account of entrainment effects between the superfluids. We find that the inertial modes are split 
into two families, ordinary fluid inertial modes (i°-mode) and superfluid inertial modes (i°-mode). The two superfluids in 
the core counter-move for the i*-modes. For the i°-modes, the dimensionless frequency acq is only weakly dependent on the 
entrainment parameter -q. For the i^-modes, on the other hand, |ko| almost linearly increases as the entrainment parameter 7] 
is increased. Because of the different dependences of kq on the parameter -q, avoided crossings as functions of rj are common 
between i°- and i°-modes. We also find that some of the i^-modes that are unstable against the gravitational radiation reaction 
at 77 = become stable invalidating the instability criterion when 77 becomes larger than some critical rjcrit, the value of which 



10 S. Yoshida and U. Lee 



depends on the mode. The radiation driven instability associated with the current muhipole radiation is quite weak for the 
inertial modes. Since the mutual friction damping in the superfluid core is quite strong for a typical value of ri compared with 
other damping mechanisms considered in this paper, we conclude that the inertial modes are easily damped by the mutual 
friction. 

Recently, Arras et al. (2002) considered non-linear couplings between an unstable r-mode and stable inertial modes of 
rotating normal fluid neutron stars as a mechanism that limits the amplitude growth of the r-mode, and discussed in the WKB 
limit the modal properties of the inertial modes, which are important to determine the saturation amplitude of the r-mode. 
In their analysis, however, no effects of buoyancy and possible core superfluidity on the inertial modes are correctly included, 
which substantially affect the modal properties such as the frequency spectra and damping rates (see Yoshida & Lee 2000b). 
If we consider cold neutron stars in a neutron star binary, inertial modes in the superfluid core of the stars will be excited by 
the tidal interaction with the companion, and the mutual friction dissipation produced by the excited inertial modes will play 
an important role in the energy and angular momentum transfer between the orbital motion and the neutron stars, since no 
p-modes are propagating in the superfluid core (Lee 1995, Andersson & Comer 2001). In fact, the damping timescale of the 
inertial modes propagating in the superfluid core may be given by r T^ip/Q, sec, and it is for Q. ~ 0.1 of order of 0.1 to 
1 sec, which is extremely short compared with the damping timescales of g-modes due to viscous processes and emission of 
gravitational waves (see, e.g., McDermott, Van Horn, & Hansen 1988; Reisenegger & Goldreich 1992). Considering that the 
orbital period decrease dlnP/dt due to gravitational radiation reactions is proportional to Z)~'* (e.g., Shapiro & Teukolsky 
1983) but that due to the tidal interaction to (e.g., Savonije & Papaloizou 1983) where D is the binary separation, the 
tidal interaction in terms of the inertial modes in the superfluid core becomes important as the binary separation decreases. 

Acknowledgments: S.Y. acknowledges financial support from Fundagao para a Ciencia e a Tecnologia (FCT) through project 
SAPIENS 36280/99. 



APPENDIX A: OSCILLATION EQUATIONS 

For the oscillation equations in the superfluid core, we employ vectors s", , t", f , , and h'' , whose components are given 
by 

s-,=Sr, sl = Sf, tl^TP, tl = Tf„ K = HT, hl = Hf, (Al) 

where I = |7n| + 2(fc — 1), /' = / -f 1 for even modes and / = |m| + 2fc — 1, = / — 1 for odd modes, where k — 1, 2, 3, • • •. 
Using these vectors, the perturbed continuity equations Q and are written as 

^ds^^^/dlnp^^gX (A3) 

dr \ dmr J 

The independent components of equations 1)26^ and H27|l are reduced to 

Lor^ - Mir^ - muAg'^r^ + {(2 + mu)l - 2muAo^}h" - (1 + 2mMo ')s^ 
dr dr dr 

+uAo^C(Ai - 21)t" + ^ {r—ih" - h") + lih" - h") - (s" - s")\ (A4) 
pn I dr J 

+4 (^) ih^-hn^o, 



dr y Pn 

Lor— - Mir— - muAg^— + {(2 + niu)l - 2niuAo^}h'' - (1 + 2miyAo^)s'' 
dr dr dr 

+i'Ao'CiAi^21)t'' + ^(l + ^^ (r-^(/i"-h'') + 2(/i"-/i'')-(s"-s'')]. (A5) 
Pp \ rrip J I dr J 

dr \ Pp \ rrip J I 

Lii" - Mo/i" + uKs'^ + £ll{e ~r)=0, (A6) 

Pn 



-1 



Lit" - Muh" + vKs" + ^ ( 1 + ) (*" - t'') = , (A7) 
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where v = 20. /a = 2/ no and 1 stands for a unit matrix. Here, the quantities, Lo, L\, C, K, Ao, Ai, Mo, and Mi are matrices, 
given as follows: 
For even modes, 



{Lo)i.i — 1 — 



l{l + l) ' 

-a + 2)J™i, 

rm 



{Li)i,i = 1- 



/ + 1 ' 

{Ao)^,^ = l{l + 1) 
I 



(/+ l)(Z + 2) ' 

= (/ + l)J,72, 

rm 
•Jl+2 



{Mo)i,i = !^ 



(Ml 



l + l 

I + 2 



■Jl+1 1 

Jm 
I + l 1 



1 + 2' 

(Ai),,, = (/ + !)(; + 2), 
/ + 3 



l + l 

where I — \m\ + 2i — 2 for i 

1/2 



(Mi)h 



/ + 2 

^Z + 2 
and 



^1+2 1 



J, 



(/ + m)(/ - m) 
(2Z-1)(2Z + 1) 

For odd modes, 

{Lo)i,i ~ 1 



(A8) 



/(/ + !) ' 

{C)^,, l)jr , (C),,, + l 



(Li)«,> = 1 



/(/-I)' 
(« + 2)J,7i, 



(Ao)i,i 



Jilt 
I 



(K) 



J, 



l{l- 



(Mo)m = jr , 



l + l' 

(Ai),,, 

(Afo)i+i,i = V 



l + l 



{Mi)i,i = u 



I 



■Ji 



(Mi)i,i+i = "^'+1 ' 



where / = \m\ +21 — 1 for i = 1, 2, 3, . . .. 

The oscillation equations in the superfluid core which we use are given as a set of linear ordinary differential equations 
for the variables and with X = n,p, which are obtained by eliminating the vectors in equations (A2)-(A5) using 
the algebraic equations (A6) and (A7). The oscillation equations in the normal fluid envelope can be obtained by substituting 
_ _ gO ^ _^p_^o — — p jjj^Q -(-jjg oscillation equations in the superfluid core, where s° and t° means the 

fluid velocity fields in the normal fluid envelope. 
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Table 1. kq for the i°- and i°-modes for the rj = model 



lo - \in.\ 


i°-mode 


i*-mode 


i-mode 


1 






0.6667 


2 


1.0978 


1.1187 


1.1000 




-0.5685 


-0.5115 


-0.5566 


3 


1.3564 


1.3865 


1.3578 




0.5180 


0.5078 


0.5173 




-1.0293 


-0.9729 


-1.0259 


4 


1.5191 


1.5522 


1.5196 




0.8639 


0.8611 


0.8630 




-0.2738 


-0.2421 


-0.2753 




-1.2738 


-1.2242 


-1.2729 



Table 2. kq and r?'s for the i°- and i^-modes of m = 2 



Modes 




n 


kq 




GJ,3^ ' 


t9, r . (s) 




rO(s) 




i° 


2 





1.0978 




-4.1 X 10* 




-2.4 X 10" 


2.5 X 


10** 


oo 








1.1187 




— 2.5 X 10 




— 1.2 X 10 


1.2 X 


10" 


oo 




3 




1.3564 


-2.2 X 10^ 




-1.7 X 10^ 




3.3 X 


10" 


oo 








1.3865 


-3.7 X 10^ 




-3.5 X 10^° 




5.2 X 


10^ 


oo 


i" 






0.5180 


-4.5 X 10^ 




-2.1 X 10« 




2.6 X 


106 


oo 


i" 






0.5078 


-5.9 X lO'' 




-1.8 X lO'' 




5.7 X 


10^ 


oo 


i° 


4 




1.5191 




-2.1 X 10" 




-2.4 X 10" 


3.8 X 


10*5 


oo 








1.5522 




-5.9 X 10^2 




-1.5 X 10"* 


2.2 X 


10^ 


oo 


i° 






0.8639 




-1.8 X 10^ 




-1.5 X 10"' 


9.9 X 


10*5 


oo 


i' 






0.8611 




-2.0 X 10* 




-5.0 X 10" 


2.3 X 


10^* 


oo 


i° 


2 


0.02 


1.0912 




-1.5 X 10'^ 




-2.1 X 101^5 


6.5 X 


10^ 


9.1 X 10-1 








1.4333 




-1.8 X 10^ 




-9.5 X 10" 


6.1 X 


10^ 


1.8 X 10-1 


i° 


3 




1.3591 


-2.1 X 10^ 




-1.7 X 109 




2.8 X 


10^ 


3.8 X lOi 


i" 






1.7829 


-1.2 X 10^3 




-1.4 X 10" 




1.6 X 


10^ 


1.9 X 10-1 


i° 






0.5171 


-4.1 X 10*^ 




-1.9 X 10« 




8.3 X 


10" 


1.4 X lOi 








0.6464 


-3.2 X 10** 




-1.8 X 10" 




3.0 X 


lO'' 


2.4 X 10-1 


i° 


4 




1.5239 




-1.4 X 10i3 




-2.7 X 10" 


1.6 X 


10^ 


3.3 X 10° 


i" 






1.9832 




-8.9 X 10^2 




-2.3 X 10'^* 


1.8 X 


lO-""' 


1.3 X 10-1 


i° 






0.8640 




-1.7 X 10^ 




-1.5 X 10"' 


2.3 X 


107 


1.4 X 102 








1.1063 




-1.2 X 10^ 




-8.9 X 10" 


1.8 X 


lO"' 


2.3 X 10-1 




2 


0.04 


1.1015 




-3.5 X 10* 




-4.9 X 10" 


2.6 X 




2.2 X lOi 








1.7632 




-4.1 X 10^2 




-8.4 X lO^s 


4.7 X 


10^ 


4.9 X 10-2 


i° 


3 




1.3605 


-2.2 X 10^ 




-1.7 X 10^ 




3.0 X 


107 


7.0 X 10° 


i' 






2.1661 


1.2 X lO^"* 




5.9 X lOi* 




1.3 X 


10-' 


4.5 X 10-1 


i° 






0.5172 


-4.4 X lO'* 




-1.9 X 10'5 




4.1 X 


10« 


4.2 X 10° 








0.7831 


-3.0 X 10^ 




-2.4 X lO^" 




1.9 X 


10^ 


6.7 X 10-2 


i° 


4 




1.5188 




-9.8 X 10^° 




-2.4 X 10" 


4.0 X 


10** 


3.3 X 10° 








2.4083 




2.2 X 10^3 




8.9 X lO^'' 


6.3 X 


10^ 


3.7 X 10-2 








0.8642 




-1.9 X 10'^ 




-1.5 X 10"> 


2.2 X 


10^ 


3.9 X lOi 








1.3476 




-1.2 X 10^° 




-2.2 X lO^s 


1.1 X 


10^ 


5.3 X 10-2 
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Figure 6. Same as Figure 5 but for even parity inertial modes having negative ko- 



Inertial modes of neutron stars with the superfluid core 17 




Figure 8. Same as Figure 7 but for odd parity inertial modes having negative kq- 
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